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Abstract 

The spectrum of anyons confined in harmonic oscillator potential 
shows both linear and nonlinear dependence on the statistical param- 
eter. While the existence of exact linear solutions have been shown 
analytically, the nonlinear dependence has been arrived at by numeri- 
cal and/or perturbative methods. We develop a method which shows 
the possibility of nonlinearly interpolating spectrum. To be specific 
we analyse the eigenvalue equation in various asymptotic regions for 
the three anyon problem. 



1 Introduction 



Non-relativistic quantum mechanics in two spatial dimensions admits possi- 
bility of fractional statistics |], |2|. Existence of fractional statistics is also 
intimately tied to having multivalued wave functions which occur naturally in 
the quantum mechanics on multiply connected spaces |J. The classification 
of multivalued wave functions is provided by one dimensional representations 
of the fundamental group of the multiply connected configuration space. Al- 
though there are several known examples of a kinematic classification Q (i.e., 
representations of the fundamental group), very little is known about the dy- 
namics on such spaces. A system of "Anyons"- particles in two dimensions 
obeying fractional statistics- is a good and physically relevent example of 
quantum mechanics on multiply connected spaces. It is physically relevent 
because anyons have been proposed as candidates for explaining fractional 
quantum Hall effect []|] and are likely candidates for high T c superconduc- 
tivity 0. It is a good example because at least some exact solutions to 
the energy eigenvalue problem for N-anyons in some external confining po- 
tential are known - an exception in many-body quantum mechanics with 
non-seperable Hamiltonians. 

Concentrating on anyons in a harmonic oscillator potential, one notices 
two important features: All the exactly known energy eigenvalues are linear 
in the statistical parameter a for arbitrary number of anyons [0]. Whereas in 
the case of three and four anyons, numerical studies exhibit some eigenvalues 
with non-linear dependence on a as well as level crossings ||. Meanfield 
studies for large number of anyons also exhibit the same feature ||. There 
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is no analytical approach exhibiting such nonlinearly interpolating energy 
eigenvalues- raising perhaps a doubt that the nonlinearity may be an artifact 
of numerical computations which are necessarily limited by computing power. 
In the present work we develop an analytical approach which does show the 
possibility of nonlinearly interpolating spectrum. We believe that the method 
is general enough to have validity even beyond the domain of the problem 
on hand. 

Our approach is based on analysing the eigenvalue equation (differen- 
tial equation) in various asymptotic regions. Noting that eigenvalues them- 
selves are independent of where the differential equation is solved, we expect 
the asymptotic form of the differential equation to determine the spectrum. 
Presently we have analysed in detail the case of three anyons in an oscillator 
potential but we believe the method can be suitably adapted to many anyons 
and possibly to some other many body systems. 

The three anyon problem in the asymptotic regions we consider reduces 
to solving a Confluent Hypergeometric Equation (CHE) in the region R < 
x < oo. This equation has a solution regular at x = and a solution which 
is irregular at x — 0. The regular solution leads to linearly interpolating 
energies while the irregular one (R is strictly greater than zero) leads to 
nonlinearly interpolating solutions. We give two first order equations which 
specify the spectrum exactly. 

The paper is organised as follows: In Sec. 2 we discuss the preliminaries, 
specify our formulation and set the notation. This part is well known and 
has been included to make the paper self contained. In Sec. 3, we give a quick 
review of known exact solutions. We point out how nonlinearity can come 
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about and note additional asymptotic regions available to extract information 
about the spectrum. After specifying the asymptotics we derive the equations 
determining the spectrum. We note how the numerically known low lying 
eigenvalues fit into our equations. Full analysis of these system of equations 
is under investigation. Section 5 contains summary and a brief discussion. 

2 Preliminaries 

In the following by "anyons" we mean a quantum mechanical system of N 
particles moving in two dimensions with wave functions which have a stip- 
ulated multi-valuedness. To make this explicit let us denote a generic wave 
function as ip(fi, fjy). Let denote the operation of taking the ith 
particle coordinate around the jth coordinate along a closed path 7. The 
path 7 does not enclose any other particle coordinate and is taken in an 
anticlockwise sense, say. Then under such an operation ip acquires a phase 
namely, 

Pijjipin, ...r N ) = exp(i2Tra)ip(f 1 , ...f N ). (1) 

If a path 7 encloses other particle coordinates as well then such a path can 
be broken into a set of closed paths each of which encloses exactly one particle 
coordinate. Applying the stipulation above, one can compute the total phase, 
for such a path. If the sense of the path is reversed then a — > —a. Clearly 
the phase acquired depends only on the homotopy class of the path (i.e., is 
same for two paths 7 and 7' if 7 and 7' can be continuously deformed into 
each other). 
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Let us introduce the complex notation for particle coordinates: Zj = 
Xj + ii/j, Zj = Xj — iijj coordinates. Clearly Zy , where Z^ — Z i — Zj, has 
the property that if Zj is taken around Zi, Zf-- changes by exp(i2na). This 
allows us to write any generic wave function satisfying eqn(l) as, 



ip(Zi,Zi) 



i<j 



4>(Zi,Zi), 



(2) 



with the bracketed expression being a phase and now ip(Zi, Zi) is a single 
valued function. 
Clearly, 



v k ^z ^ ,z■ i ) = X{&) a,2 

i<j 

which can be rewritten as, 



V k ^(Z h Zi) + V fc /n(n(^) a/2 )^(^, ^ 



(3) 



v fc ^,z l ) = n(yr /2 

i<j ^ij 



V k *P(Zi, Zi) + ia J2 7 * r *2 ^( z ii z i) 

j+k I r kj I 



(4) 



Since ip is single valued, the r.h.s of the above equation has exactly the same 
multivaluedness as the l.h.s. In otherwords we have, 



where 



and 



V4n(ft a/2 ^ Zi)] = Dk ^ 

i<j i<j 



D k ij = [V fe + A k ]4> 



(5) 



j+k I r kj I 

This allows us to write any higher order differential operators on ip in 
terms of corresponding covariant differential operators on the singlevalued 



wave function tp. In particular a Hamiltonian operator, typically — J2i Vf + V 
can be written similarly. An eigenvalue equation written in terms of ip can 
then be recast as a corresponding equation in terms of ip. 

Although both formulations are equivalent, dealing with operators on 
multivalued wave functions is much less transparent than dealing with oper- 
ators on singlevalued wave functions. Naive commutation rules, symmetries 
that one would expect by looking at an operator on single valued functions 
are not at all true in general for the "same" differential operators acting on 
multivalued wave function. 

Considering eigenvalue problem in terms of ip has other advantages too. 
Since all the subtleties of multivaluedness are equivalently transcribed in 
terms of additional "interaction" terms (the so called statistical interactions), 
the eigenvalue problem is amenable to approximations. One is also on firmer 
ground in doing usual algebraic manipulations with operators. With these 
in mind we will work with singlevalued wave functions with "statistical in- 
teractions" . 

As a first step one would like to understand the system of "free anyons" . 
However, the statistical interaction falls off So 
one is not sure whether the Hamiltonian with only statistical interactionss 
has only discrete eigenvalues. One can put the system in a box to ensure 
discrete eigenvalues but then one needs suitable boundary condition. An 
oscillator potential ensures discrete spectrum without introducing a finite 
size. One could take some other confining potential but in the limit a — > 
one should know the spectrum. One then has hope of doing at least the 
perturbative analysis. Since the statistical interaction depends only on rel- 
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ative separations, the Centre of Mass(CM) dynamics should play a trivial 
role and oscillator potential also allows a separation of CM and relative co- 
ordinate dynamics. The oscillator potential problem can also be mapped on 
to a problem of anyons in a real, constant external magnetic field along the 



Z-axis ||10|| . Bearing these facts in mind, we choose the oscillator potential 
without further justification. In order to derive the thermodynamic proper- 
ties of a system of anyons there exist well defined methods of eliminating the 
dependence on the oscillator frequency |T|. 

The Hamiltonian we consider is- after carrying out the usual scaling of 
variables (h = c = 1) - interms of dimensionless quantities, 

1 N 1 N N £■■ a 2 N r-f-u 

H = hu[W^ + l^-a £ ? + yES ( 6 ) 

z i=l * i=l j>i=i'ij ^ ijtj,k ' ij' ik 

where 

kj = (ri -rj)x (Pi-Pj). 



and all distances have been expressed in units of 1/ y/muj. Notice that the sta- 
tistical interaction is independent of the centre of mass. This is the operator 
we analyse subsequently. 



3 Asymptotic Analysis 

First, we briefly discuss the class of exact solutions already known and then 
move over to the asymptotic analysis. For discussing the known class of exact 
solutions it is convenient to use the complex coordinates Zj, Zi in terms of 
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which the Hamiltonian takes the form, 

H = ~2dA + \z iZl - a £ A ~ ^) + % E (7) 
where <9j = d/dZf, = d; t — dj, etc., and the eigenvalue equation is, 

Hi/>(Z i ,Z i ) = E^(Z i ,Z i ). (8) 
The conserved angular momentum J, is given by, 

J = ZA - ZA (9) 

For the sake of completeness we recapitulate the standard asymptotic anal- 
ysis, i.e,Zi = \ Zi , A — > oo : 

(i) As Zj — > oo, Zy — > oo, the oscillator potential term dominates over the 
a dependent terms reducing H to the Hamiltonian of N two dimensional 
oscillators. All the eigenf unctions of the oscillator have the Gaussian 
suppression factor. Thus we put, 

rl>(Z i ,Z i )=exp(-^Z i Z i )MZ i ,Z i ). (10) 

Substitution of this form transforms the original eigenvalue equation to an 
equation for ipi, viz., 

l-2dA+ZA+ZA-aY,&-^) + ^- E yV^ 1 = ( E ~ N ^ ( n ) 

i<j ^ij ^ij A ijtj,k ^v^ik 

The N on the r.h.s is the zero point energy for N two dimensional 
oscillators while on the l.h.s the oscillator potential is traded for the 
"scaling operator", ZA + ZA. 
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(ii) The other asymptotic region is | Zij | — > for any pair of particles i and 
j, say. Since the a dependent terms are singular we expect the wave 
function to vanish sufficiently fast in this limit so that the eigenvalue 
equation is well defined. 

Assuming a Taylor series expansion in as | Z^ | — > , we can write, 



for A > and we have taken the pair (12) for definiteness. Substituting this 
in the eigenvalue equation and keeping terms to the leading order in | Z 12 \ 
gives, 



or A = ±<x However A > 0, thus A = a. Since this must happen for each 
pair of particles seperately and independently we can take out a factor 
| X | a , where 




(12) 



a > 



A 2 - a 2 = 



x = l[z tr 



(13) 



i<j 

For future reference let us also note that if 



'2, 



(14) 



then the eigenvalue equation transforms to 



[-2^ + ^ + ^ ! -(a + 7)E| i + («-^)E^ + 



i<j ^ij i<j ^ij 





N(N-l) fa 



Taking (3 = 7 = a, we have, 



fa=\X\ a fa, 



(15) 
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and 

_ _ _ fi- ■ ry 

[-2dA + ZA + ZA - 2aJ2 it } ^2 = [E - N - -N(N - 1) ] ^2 (16) 
while if j3 — 7 = —a, we have, 

^=|X|- a ^, (17) 

and 

— — — &■ ■ Ql 

l-2d t d t + ZA + Z l d t + 2aJ2^ L ]^ 2 = [E-N + -N(N - 1) ] ^ (18) 

i<jr »J 

Clearly in the second case ip' 2 will have to vanish at least as fast as | X \ 2a . 
Also, if Jipi = , then for general (3, 7 we have 

^■02 = bi 2 2~ ^ = ' 

For the two cases noted above, /3 = 7 = ±a, j 2 = j\. Incidentally if (3 — 
a, 7 = —a choice is made, then 

^ = (§r /2 ^' (19) 

and 1P2 satisfies the simple equation, 

[-20A + ZA + ZAWi =[E- N]M. (20) 

That is we recover the so called "anyon gauge" where the wave functions are 
multivalued. If we want to retain the singlevaluedness of the wave function 
then our choices are limited to 7 = f3 + integer. The three-body interaction 
term vanishes if 7 = —aor-j3 = a. Allowing a non-zero integer only changes 
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the angular momentum of the wave function with a proportionate change in 
the energy. We will, to begin with retain the angular momentum also and 
thus have (5 = 7 = ±a as the two choices. 

At this stage we can summarise the known exact solution in a convenient 
way @. With the choice f3 — 7 = a, ipi —\ X \ a ip 2 , we have two basics types 
of exact solutions: 

with the energy eigenvalues given by, 

r, »r N(N-l) 

E = N - ? + a— - 

J 2 

and 

^ = <Mt = El^l 2 ); i = 0- ( 23 ) 

i 

with ip2{t) a polynomial of degree m in t. The corresponding energy 
eigenvalues are given by, 

E = N + 2m + a iV ( iV ~ 1 ) i (24) 

The second solution is necessarily bosonic since t is symmetric, where as 
the first solution needs explicit symmetrization and antisymmetrization 
of the wavefunction in terms of Zi to obtain the bosonic and fermionic 
wavef unctions. Since this is always possible, the degeneracy of the first 
type solution is exactly the same for both bosonic and fermionic type 
solutions for any given angular momentum j (< 0); 
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(21) 
(22) 



Obtaining the first solution is trivial since only the "scaling operator" 
contributes and —j in the energy eigenvalue is nothing but the de- 
gree of scaling. The second solution can be obtained most easily by 
transforming the equation in the form of a confluent hypergeometric 
equation: 

[t^ 2 +(b-t)f t -a]Mt)=0, (25) 

where 

Irrn »r N(N~l), 

a = — [E — N — J -} 

2 l 2 J 

b = N + a — ^— 

The admissible solutions are the regular Confluent Hypergeometric 
Functions(CHF), M(a,b,t) . That a = —m follows from the demand 
of normalizability of the wave function. The corresponding eigenvalues 
follow immediately. 

We may also take a combination (product) of the solutions the types 
discussed above, Eqns(21) and (23), to get further j < solutions. 
These are : 

fo = f(Zi)g(t); E mJ — 2m — j + N + a N ^ N ~ ^ . (26) 

• With the choice f3 — 7 = —a and ipi =| X \~ a ^ we have only 
the analouge of the first type of solutions i.e. ip2 = ^{Zj). ip2 has 
to be of the form X d f(Zi) with d > 2a so that ipi vanishes if the 
coordinates of any two particles coincide. This means that j has to be 
large enough(j > (N — 1)(N — 2)/2). The corresponding eigenvalues 
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are given by 



If however the condition, j > (N — l)(iV — 2)/2, is not satisfied then 
the wave function ipi remains regular only for some values of a (0 < 
a < 2j/N(N — 1)) but not for all < a < 1 which gives rise to the so 
called noninterpolating solutions which have also been discussed in the 
literature Jl3|. To the best of our knowledge these are all the exactly 
known solutions for the N anyon problem. 

All these solutions for the energy eigenvalues have a linear dependence on 
a while the corresponding eigenfunctions are finite order polynomials apart 
from the | X \ a and the Gaussian factor. In fact a simple scaling arguement 
shows that if an exact eigenfunction is a polynomial (i.e. has a finite degree), 
then the corresponding eigenvalue must be linear in a. 

For, if ip2 is a polynomial in Zi,Zi, with highest total degree d then 
for Z i} Zi — > oo the polynomial becomes a monomial and only the scaling 
operator term dominates giving, 

E = N + d + a N( - N ~ X \ (28) 

Now if | Zij |— > every eigenfunction vanishes as | | a or faster, i.e., after 
| X \ a is taken out the remaining function has only integer powers of Z^, Zi. 
Thus d has to be an integer and hence E has a specific linear dependence on 
a:. 

However it is well known that solutions which have nonlinear dependence 
exist. This has been shown numerically ||, through perturbative analysis 
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for three anyon problem and using meanfield methods for large number 
anyons ||. This point is also our focus from now on. The scaling arguement 
will then imply that d has to be a nonlinear function of a. How is this 
achieved? 

Suppose ip2 admits a power series representation with infinite radius of 
convergence but the series does not truncate then the scaling arguement fails. 
However the analysis of such a series solution always seems to lead to ex- 
ponentially divergent behaviour making the solution non-normalizable, ie., 
the power series has to truncate. But then only linear soutions are possible. 
The other possibility then is that the power series has a finite radius of con- 
vergence, i.e., a "scale" R has to enter if nonlinearity is to be possible. But 
this means that ip2 has two different representations as the scale parameter 
A — > and A — > oo, each being valid for A < \ max and A > A m f„ respectively. 
For A > Xmin the scaling arguement can still work but now the series need 
not have only integer powers, i.e., d can be a nonlinear function of a. We 
conclude then that if a nonlinearly interpolating eigenvalue is possible at all, 
the corresponding eigenfunction must have two different series representa- 
tions for ^(AZj, XZi) as A —>■ and A — > oo. One then has to try to match 
the two series suitably. 

Exploring this possibility looks extremely complicated for general N. So 
we now specialise to N=3 and go over to relative coordinate dynamics to re- 
duce the number of independent variables in the differential equation. These 
are given interms of the Jacobi coordinates, 

Zi — Zo , , 
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P 2 = j= 3 (30) 

and similarly their complex conjugates. 

After separating out the CM coordinate and using the notation pi = 
TiexpiiOi) the relevant hamiltonion becomes, 

3a r r| - T \ e ^-^) d id. 



„ / <9 id,, 

+2r 2 — 

<9r 2 r 2 <90 2 

The full wave function is of course, 

^ =| X | a exp(-(l/2)^r2)0(r 1 ,r 2 ,^ 1 ,^ 2 )$ CM ( J R,e). (31) 

i 

So far two asymptotic regions have been used to simplify the original 
eigenvalue equation: 

(i) r 1 ,r 2 — > oo such that r±/r 2 is fixed and finite. 

(ii) ri,r 2 — > such that ri/r 2 is fixed and finite. 

and we arrived at the possibility of two different series for some wave func- 
tions. Clearly there are two more regions possible, namely, 

(iii) r\ — > such that r 2 is fixed or r 2 — > oo such that r\ is 
fixed, i.e.,ri/r 2 — > 0. 

(iv) r 2 — > such that r 1 is fixed or ri — > oo such that r 2 is 
fixed, i.e.,r 2 /rx — > 0. 
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Notice that these types of regions become possible only for iV > 3, when 
there are at least two relative coordinates. For two anyon problem there is 
a single relative coordinate and asymptotics for <fi can be considered only as 
r — > oo or r — > (regions (i) and (ii)). Since two anyons have no possibility 
of any finite R or scale parameter entering one has to use Taylor series, and 
hence the wave functions are necessarily polynomials. This leads to a linearly 
interpolating spectrum. This is true provided either that a dependence is 
assumed to be smooth or is C°° function. 

To focus, let us consider (iii). As r 1 /r 2 — * the Hamiltonian becomes 
separable and Ji = —i-^- commute with the corresponding separated Hamil- 
tonian. This is true under the assumption that the derivatives in the a 
dependent terms in the H re i remain finite as ri/r 2 — > 0. 



That Pi — a is obvious because by definition r\ involves particles 1 and 2. 
P 2 = 2a follows from the fact that r 2 involves all the three particles, and 
in the asymptotic region particle 3 can only wind around particles 1 and 2 
together, but not around any one of them seperately. 
Because of this as r 1 /r 2 — > 0, 




where 



Pi = a; p 2 — 2a 



and 



H = H 1 + if 2 . 




(33) 
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where, 

<t>i{n, 0%) = r] i (r i )exp(in i 9 i ); n 1 +n 2 = j. 
Now consider ri/r 2 — > as ri — > and r 2 finite and fixed. Let, 

</>(ri,r 2 A,0 2 ) — > r*exp(in 1 9 1 )(f) 2 (r 2 ,e 2 ). (34) 

Since H is seperable, we may first compute ffi acting on 0. We have, 

Htf = r^p(mA)0 2 (r 2 A)[A - \ ^ " + 2 ° (Wl + A)) ]. (35) 

Because ri — > limit has to be finite (leading power is already taken out) we 
get, 

A 2 — n\ + 2a{n\ + A) = =>- A = — a± | ni — ok | 

Now in the same limit , we have | X \ a -^ rfr 2 a and therefore the total power 
of T\ becomes ± ] n\ — a |. Since has to vanish as r\ — > (r 2 fixed) we 
have to choose the positive sign. 

This also matches with the two anyon wave function in relative coordi- 
nates which vanishes as rl n_a L Thus choosing the positive sign, 

A = — a+ | ni — a | 

we have 

{H 1 + H 2 )(j> = [(H 2 + \)<t) 2 \r$exp(in 1 9 1 ) (36) 
= (E-2- 3a)(j) 2 r^exp(in 1 9 1 ) (37) 

Thus (j) 2 satisfies the eigenvalue equation 

# 2 </> 2 = (E-2 -3a- \)<f) 2 = [E-2- 2a- \ m - a \)<j> 2 = S</) 2 . (38) 
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H 2 is once again the "two anyon Hamiltonian" . There is one crucial differ- 
ence; since our equations are valid only asymptotically r\/r 2 — > we have 
to maintain r 2 nonzero. Let us say that r 2 > R , where Rq is some new 
parameter. 

With the definitions 

H = -[-2a+\n 2 -2a\] (39) 
b = 1+ | n 2 - 2a | (40) 
a = --[£ + 2a- | n 2 -2a \ ] (41) 

= --[E-2- | ni -a \ - \ n 2 -2a \ } (42) 
x = r\ (43) 
R = Rl (44) 



>2 



r 2 ,9 2 ) = x^exp(in 2 9 2 )u(x) (45) 



the eigenvalue equation for 2 becomes the CFE for u(x): 

d 2 d 
[x— + (b-x)—-a]u(x) = 0, (46) 
dx l dx 

R < x < oo. 
The general solution is given by, 

u(x) = C x M{a,b,x) +C 2 U(a,b,x), (47) 



where M and U are in standard notation [TJ| and C\ and C 2 are constants. 



For C\ we have two possibilities, namely it is zero or non zero. If it 
is non zero then normalizability as x —>■ oo requires a to be minus a non 
negative integer. This is precisely the case when the U solution ceases to be 
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an independent solution. In other words, if C\ is non zero then we really 
have only one solution which is a finite order polynomial. This solution is 
also well behaved as x goes to zero and therefore we may take R to be zero. 
The energy is read off by putting the value of a in the definitions given above. 
We see that the energy eigenvalue in this case is linear in a. The known exact 
solutions are covered by this case. 

If C\ is zero then C2 is non zero and we have only the U solution. U is 
well behaved for large x and thus we do not get any condition on a from the 
demand of normalizability. However U is divergent as x — > and we have to 
keep R strictly positive. 

The only way E can get quantised now is by putting a boundary condition 
at x — R. Hermiticity analysis implies U \r= or dU/dx \r— 0. But 

— — = —all (a + 1, b + 1, x) 
dx 

Because a is not a non positive integer (else we are back to the previous case 
) essentially we get the conditions 

U(a,b,x) = 

or 

U{a + 1,6 + 1, x) = 0. 



Now U(a,b,x) has real positive zeros only if [[14 

a < and 1 + a — b < 0. 
In the present case b = 1+ | n<i — 2a \ and 

1 + a — b = a — \ n 2 — 2a \ 



Zeros of U exist only if, 

a < and a < | n 2 — 2a \ (if U = 0) 
or 

a + l<0 and a < \ n 2 - 2a \ (iflf' = 0), 

that is for a < or a < -1 respectively. In the limit a^Owe know that all 
solutions are finite polynomials. In effect then R — > must hold as a — > 0, 
i.e., a must be a negative integer. The Bosonic ground state has energy 2hu. 
To admit this possibility we must allow for a to be zero. We will therefore 
choose U (a, b, R) = as our boundary condition. Since for a < (a not an 
integer) U has in general more than one zero we will stipulate R to be the 
smallest zero of U (say). This then determines uniquely a and hence E as a 
function of R which is still an unknown parameter. 

It may seem that the above boundary condition is very restrictive. It 
is not. For suppose we knew E = E(a), then U has zeroes and hence 
the smallest zero. Therefore for all a, we can always find an R such that, 
U(a(a),b, R) = 0. But now clearly R depends on a and the asymptotic 
quantum numbers ni,n 2 . Thus a(a 1 n\ 1 n2) determines R = R(a,ni,n 2 ). 
Conversely, if we could determine R = R(a,ni,n 2 ) by some arguement then 
U = will determine E^a, n 1? n 2 ). 

The boundary condition can thus be understood as stipulating that there 
exists functions a(a),R(a) and b(a) which is explicitly known, such that 
U(a,b,R) = 0, i.e., in the space parametrised by variables E, b and R with 
2<_E<oo, 6=1+ | n 2 — 2a ], < R < oo there exist curves along which 
U=0. Determining these curves is equivalent to getting the spectrum. The 

19 



differential statement of U = along a curve is of course = which can 
be written as: 



dU da 



dU db dU dR 
db da dR da 



(48) 



da da 



+ 







Meanwhile we have tested the U = criterian in the following way. We 
choose a particular state at the fermionic or the bosonic end. Use the nu- 
merically determined eigenvalue to get a(a) and determine the zero of the U 
function numerically. Thus we get a R(a). We fit a form to thus determined 
R and the compute E at arbitrary values of a. This still has predictive power 
in that once we postulate a form for R(a) we can determine its parameters 
for small values of a where perturbation theory can be used. We then have 
a prediction for all values of a in the range [0,1]. In particular we have car- 
ried out this exercise for the state which interpolates to the fermionic ground 
state and the corresponding bosonic state with j=2 which is supersymmetric 
M to the fermionic ground state. The results are given in the figure 1 for 
the energy and figure 2 for the corresponding R . 

In the end we briefly mention another method of analysing solutions in 
the three anyon problem using the Fourier expansion of the wavefunction 
in eq.(31). Note that 4> may be written as, 



for a given j, where j = n\ + ni, n = n\ — n^- We can now use this 
representation to look at the solutions for a given j. The eigenvalue equation 
now becomes H re i<j)j = S<f)\, where H re i is the three-anyon Hamiltonian in 



n 



(49) 
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relative coordinates and S = E — 2 — 3a as before. Substituting for the 
eigenvalue equation takes the form, 

or. 92 9 (j+n)(j+n-4a) 
3rr 2 4ri— ^ + 4(1 + a - Xi)- 

<9 2 w- x <9 (?' — n)(? — n — 8a) 
+ 4^ + 4(1 + 2.-,,)--^^- <+2£ te „ 

T4 <9 2 „/, o n 5 (j + n + 4)(j + n + 4- 12a) 

<9 2 9 f?'-n-4) 2 

+ 4x^ + 4(1-^)—- ^ - +2£1 W4 , 

where = r 2 . This equation is exact. These infinite set of coupled equa- 
tions relate Fourier modes differing by 4, i.e., if xo, Xi; X2, X3 are know then 
X4fc,X4fc+i,X4fc+2,X4fc+3 get determined in terms of Xo,---,X3 etc. However 
there is no relation among Xo> Xij X2 5 X3 themselves. In a sense these four 
functions will give four independent solutions of the eigenvalue equations. 
We can then deal with a given "tower" separately and independently and 
this is true for every given j. Let us concentrate on one tower. Now three 
distinct cases arise naturally: 

(a) Only one member of the tower is nonzero, ie., Xj,n = Xj,m8 n ,m- 

(b) Only a finite number of Xj,n s are nonzero, ie., 
Xj,n = V n > rii and V n < n 2 with n 2 < n\. 

(c) Infinitely many Xj,n$ are nonzero. This gives raise to the following cases: 

Xj, n = V n > n x 
Xj, n = V n < n 2 
Xj,n ^ V n. 
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Case (a) is simple to analyse and it reproduces the known exact solutions. 
In the exact eigenvalue equation relating Xj,n an d Xj,n+4 given above we set 
lhs to zero which gives an equation for the Xj,n- But the rhs also gives rise to 
another equation for the same Xj,n when the coupling to Xj,n-4 is taken into 
account. However the form of equations suggests that the two equations are 
seperable but must be consistently solved. The consistency conditions imme- 
diately yield the linear exact solutions which are already outlined previously. 
Cases (b) and (c) seem quite complicated and nonlinearly interpolating states 
must be in one of these cases. Consider the normalisation condition on the 
full wavefunction, 

roo _ 

II ^ f = E / d Xl dx 2 \X\ 2a \ Xj ,n\ 2 e- Xl+x ' =Y,°n (50) 

If ipj is given by case (b) then there are only a finite number of terms in the 
norm and if C ra 's are finite then ipj is normalisable. Thus the quantisation 
condition for the nonlinear states must arise by the demand that C n must 
be finite. For the case (c) even if all C n 's are finite we may still get the 
wavefunction to be infinite because the sum above may not converge and 
then the quantisation condition would be the convergence of I]„C n . 

We have not at present analysed this scenario in detail. But it does 
appear that the case (b) is ruled out in which case the exact solutions come 
from case (a) and the nonlinear solutions arise from case (c). 
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4 Summary and Discussion 



We have in this paper, analysed the problem of three anyons confined in 
an oscillator potential. The energy spectrum of three anyons generically 
contains two types of interpolations as a function of the statistical parameter 
a. The first type are the linearly interpoating eigenvalues with AE — E(a — 
1) — E(a = 0) = 3. These solutions are easily generalised to N-anyon case 
where AE = N(N — 1)/2,N > 2. The second set consists of solutions for 
eigenvalues which are in general non-linear functions of a. In this paper we 
have outlined an analytical method of obtaining such solutions. We exploit 
the known asymptotics to arrive at these solutions. Since the differential 
equation is valid in all regions it is sufficient to specify the wavefunction 
asymptotically to obtain the exact eigenvalues. 

For N=3, there are two relative coordinates whose magnitudes are given 
by 7"i,r 2 . The asymptotic regions are then defined by, 

• (a) ri,r 2 — > oo. 

• (b) r 1 ,r 2 — > 0. 

• (c) r\ — > 0, r 2 fixed; or r 2 — > oo,ri fixed. 

• (d) r 2 — > 0, T\ fixed; or r\ — > oo,r 2 fixed. 

Asymptotic region (a) immediately yields the Gaussian which controls the 
behaviour of the wavefunction at infinity and region (b) yields the Laughlin- 
Wu factor which regulates the wave function as two anyons approach each 
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other. The asymptotic region (c) (or (d)) is used to obtain the nonlinear ly 
interpolating solutions. 

There is a caveat to note: There may be exact eigenfunctions which do 
not satisfy the asymptotic equation (33). In this case our analysis can at 
best be expected to give an approximation to the exact eigenvalues. In the 
absense of a complete rigorous proof precluding this possiblity we have to 
note it as an open one. 

A few additional points are worth noting. 

i) We have chosen R to be the smallest zero of U without any justifica- 
tion. For low lying eigenvalues the U function has only one zero so there 
is no ambiguity. For higher states though more zeros are possible and their 
implications are to be explored and elaborated. 

ii) Some qualitative properties such as a tower of evenly spaced eigenval- 
ues for every non linearly interpolating eigenvalue || and "supersymmetry" 
for a class of eigenstates || have been noted in the literature cited. Whether 
and how these fit into our appraach is to be seen. 
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Figure Captions 

Figure 1: Plot of the Fermionic ground state energy and the corresponding 
Bosonic j = 2 state as a function of a obtained by using the boundary 
condition at R. The non-interacting Bosonic state is at a = and the 
non-interacting Fermionic state is at a = 1 and energy E = 4 in 
dimensionless units. The points shown along the Fermionic curve are taken 
from the numerical calculations ||. 

Figure 2: Shows the values of R as a function of a used in computing the 
energies of the corresponding states. For the Fermionic ground state 
R pa 0.405a(l — a) obtaied from numerical computation close to a = 0. 
The corresponding form for the Bosonic state is R ~ a 6 . 
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